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1.  Background 


The  roots  of  this  work  are  from  a paper  by  Richard  Simon  [1], 

"Stationary  Properties  of  a Two-Echelon  Inventory  Model  for  Low  Demand 
Items."  A subsequent  paper  by  Kruse  and  Kaplan  [2]  claimed  to  correct 
Simon's  results,  but  Kruse  [3]  points  out  that  Simon's  results  are,  never- 
theless, correct  in  spite  of  the  apparent  logical  error  in  their  derivation. 
He  also  notes  that  both  Simon's  and  Kruse  and  Kaplan's  results  reduce 
to  a simple  and  intuitive  form. 

We  take  a new  look  at  the  model  in  this  paper  with  the  intent  of  pro- 
viding a more  revealing  analysis  than  appears  in  either  [1]  or  [2].  We 
then  show  how  these  2-echelon  results  can  be  extended  to  model  an  N-echelon 
inventory  system  in  a manner  similar  to  the  way  the  METRIC  model  uses  Palm’s 
Theorem  to  accommodate  an  indefinite  number  of  echelons. 

With  apologies  to  Dr.  Simon,  we  call  this  the  Simple  Simon  model,  not 
only  because  it  is  an  apt  descriptor,  but  also  because  that  designation, 
as  we  have  learned,  is  irrepressible. 

2.  System  Description 

Simon's  original  work  was  done  for  RAND  and,  consequently,  uses  Air 
Force  terminology.  We  will  continue  that.  There  is  a lowest  echelon, 
composed  of  several  independent  bases,  which  is  the  source  of  failures 
for  a particular  item.  When  an  item  fails,  it  is  inspected  at  the  base 
and  a decision  is  made  to  repair  the  item  at  the  base,  repair  the  item 
at  the  depot,  or  forego  repair  and  dispose  of  the  carcass.  Spares  may  be 

held  at  base  or  depot.  If  the  failed  item  is  to  be  repaired  at  the  base, 

a replacement  item  is  supplied  from  the  base  spare  stock,  and  the  failed 
item  enters  the  spare  stock  after  being  repaired  itself.  Spare  stock  at 
the  base. is  issued  on  a FIFO  basis  and  all  demands  are  backordered  until 

^illfi^*  If  the  item  is  not  to  be  repaired  at  the  base,  a spare  is  still 

supplied  from  base  stock  and  a replenishment  spare  is  requested  from  the 
depot.  Demands  on  the  depot  are  also  filled  with  FIFO  priority  and  back- 
ordered until  filled.  Items  which  the  depot  repairs  go  into  its  spare 
stock.  Other  replenishment  comes  from  an  exogenous  source  with  an  infinite 
supply. 
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3.  Assumptions  and  Definitions 


For  the  reader’s  convenience,  we  repeat  the  assumptions  and  definitions 
from  [1]  and  [2]  below: 

r^  =“  the  probability  that  an  item  failed  at  base  j will  be 
repaired  at  base  j 

p « the  probability  that  a failed  carcass  that  is  not  base 
reparable  will  be  depot  reparable 

Dj(t),  t ^0  = the  demand  process  at  base  j (j  * 1,2,...J).  This  is 

assumed  to  be  a simple  Poisson  process  with  rate  . 

•Q 

Dj(t),  t ^ 0 =*  the  stochastic  process  of  units  repaired  at  base  j. 

Dj(t),  t ^ 0 = the  stochastic  process  of  units  sent  to  the  depot  from 

base  j for  repair. 
r 

Dj  (t) , t ^ 0 = the  stochastic  process  of  units  condemned  at  base  j. 

J 

D^(t)  = Z D^(t) 

° 3=1  ^ 

c c 

D^t)  = E D^(t) 

D (t)  * D^(t)  + D°(t) 
o o o 

Bj (t)  = the  number  of  backorders  at  location  j at  time  t; 
j = 0 denotes  the  depot. 

= the  deterministic  repair  time  at  facility. 

tj  =*  the  deterministic  delivery  time  from  depot  to  base  j. 

t^  * the  deterministic  procurement  lead  time  from  the  exogenous 

source  to  the  depot.  We  assume  R < t . 

o o 


Sj-l,Sj  =»  replenishment  policy  used  at  base  j. 

8 ,S  * replenishment  policy  used  at  depot, 
o o 

N(tj^,t2)  “ N(t2)  - N(tj^)  for  any  stochastic  process. 


A . Stationary  Distributions  of  Base  Backorders 


As  a consequence  of  the  inventory  policy  at  base  j we  have 


that 


Bj(t)  = X^(t)  + Yj(t)  - 


where 


and 


(t)  = amount  due  in  to  base  j at  time  t from  the  depot, 
(t)  = amount  in  repair  at  base  j at  time  t. 


Bj  (t)  < 0 indicates  stock  on  hand. 

By  virtue  of  the  assumptions,  independent  random  vari- 

ables, and  for  t > R 


^ ~X  r R 

Pr[Y^(t)  = k]  =e  J J ^ (X^rjRj)^/kI 


Now  (t)  t is  composed  of  two  mutually  exclusive  and  in- 

dependent elements  X (t)  and  X «(t)  where 

J > 1 J > A 

Xj  ^(t)  =»  amount  due  in  to  base  j of  those  demands  placed  on 


and 


the  depot  in  (t-tj,t) 


X «(t)  = amount  due  in  to  base  j of  those  demands  placed  on 

J > A 


the  depot  prior  to  t-t 


Clearly, 


Pr[X  (t)  = k]  =€ 

J 1 


j' 

-(X.) (1-r  ) (t.)  , 

J ^ ^ ((X^)(l-r^)(t^))VkI 


The  distribution  of  X „(t)  remains  to  be  found.  Note  that 

J 

X.  ^(t)  is  also  the  number  of  backorders  at  the 

! f “ 

depot  at  time  t-t.  which  are  for  base  j.  In  [I]  and  [2],  X.  «(t)  was 
J J J 

called  ^j(t)  iis  distribution  was  tediously  derived  by  analyzing  the 

composition  of  the  demand  sequences  in  the  interval  t-t^-t^  to  t-t^  as 

compared  to  the  depot  asset  position  at  time  t-t^-t^ . The  analyses  made 

a careful  distinction  between  the  demand  process  associated  with  reparable 

carcasses  and  the  process  associated  with  condemnation.  Motivating 


this  distinction  were  the  realizations  that  the  number  of  depot  backorders  at 

C D 

t-tj  were  dependent  on  only  D^(t-t^-tj ,t~tj)  + D^(t-R^-t^ , and 

in  [2]  that  D^(t-t  -t..t--R  -tj  could  affect  to  whom  those  backorders  at  t-t, 

**  0^  o j * o j j 

were  due*  However,  a seemingly  innocuous  assumption  that  the  probability 
that  a failed  carcass  that  is  not  base  reparable  will  be  depot  reparable, 
i.e.  p,  is  base  independent  allows  a dramatic  simplification  in  the  results 
as  well  as  the  analysis.  Moreover,  as  noted  in  [3],  this  assumption  ex- 
plains why  the  results  in  [1]  and  [2]  are  the  same  upon  algebraic  reduction. 

Suppose  there  is  a demand  on  the  depot  and  we  wish  to  know  the  probability 
that  its  source  was  base  j , If  the  demand  were  associated  with  a reparable 
carcass,  then  because  that  process  is  Poisson  at  each  base  the  probability  is 

X.(l-r  )p/  IX^(l-r^)p  (1-r  )/  IX,  (1-r^) 

J J k ^ ^ k 

If  the  demand  were  associated  with  a condemned  carcass,  the  probability  is 
likewise 


■ X (l-r  )(l-p)/  EX  (1-r  )(l-p)  - X.(l-r  )/  SXj^d-r^) 

So,  assuming  that  p is  base  independent  means  that  the  type  of  demand  does 
not  affect  the  probability  that  the  demand  comes  from  a given  base. 

A natural  way  to  derive  the  distribution  of  X ^(t)  is  through  its 
conditional  distribution  on  total  depot  backorders  at  t-t^ . 

That  is , 

oo 

Pr[Xj  ^it)  = k]  = E ” n]Pr [B^ (t-t^ ) *•  n] 

Now,  given  that  ^ provides  probability  information  on  the  number 

of  demands  associated  with  a reparable  carcass  and  the  number  associated 
with  a condemnation  which  are  backordered  at  • But,  by  the  argument 

above  this  information  is  irrelevant  for  determining  the  probability  a 
given  backorder  is  for  a given  base.  Consequently,  and  because  the  demand 
processes  are  Poisson,  we  have  that 


= n]  - '(JJ) 


where  P = X (1-r  )/ 


Moreover 


S 


o 


lim  Pr[B  (t)  = n]  * z 
_ o 


t— >00 


O 


+ (t-R  -t,  = a + n] 

o " 0 j * 


is  no  reduction  to  a simple  form. 

These  computations  can  be  done  quickly  by  computer  using  an  approach 


For  higher  demand  items,  these  arrays  will  overflow  and  a more  storage 
conscious  approach  is  needed. 

5.  Extension  to  N Echelons 

Suppose  we  consider  any  imbedded  two  echelon  structure  within  an 
N Echelon  inventory  system  of  arborescent  structure  (see  figure  below) 
for  which  the  above  assumptions  hold.  With  regard  to  our  assumptions, 
we  interpret  the  term  "condemnation*’  to  represent  the  situation  when  an 
item  cannot  be  repaired  at  the  top  echelon  in  the  imbedded  two  echelon 
structure  and  must  be  sent  to  the  next  echelon.  Further,  suppose  that  the 
stationary  distribution  of  backorders  is  known  for  the  upper  echelon  in 
the  imbedded  two  echelon  structure.  Then,  we  can  use  the  analysis  in  the 
previous  section  to  get  the  stationary  distribution  of  backorders  at  each 
lower  echelon  location.  This  allows  us  to  go  one  echelon  lower  in  the 
structure  whereupon  the  present  lower  echelon  locations  each  become  an  upper 
echelon  location  in  new  imbedded  two  echelon  structures. 

This  method  is  analagous  to  the  approach  used  in  METRIC  to  model 
N echelons  except  that  Simple  Simon  explicitly  uses  the  distribution  of  upper 
echelon  backorders  while  METRIC  uses  only  their  expected  value  thru  Palm’s 
Theorem. 
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THREE-ECHELON  ARBORESCENT  STRUCTURE 


FIGURE  1 
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